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Abstract 

It has been suggested recently to study the dynamics of a gravitating gluon condensate q in the 
context of a spatially flat Friedmann-Robertson-Walker universe. The expansion of the Universe 
(or, more generally, the presence of a nonvanishing Ricci curvature scalar R) perturbs the gluon 
condensate and may induce a nonanalytic term h{R, q) in the effective gravitational action. The 
aim of this article is to explore the cosmological implications of a particular nonanalytic term 
h oc rilRl^/"^ kP^^- With a quadratic approximation of the gravitating gluon-condensate vacuum 
energy density pvil) near the equilibrium value go and a small coupling constant r/ of the modified- 
gravity term h, an "accelerating universe" is obtained which resembles the present Universe, both 
qualitatively and quantitatively. The unknown component X of this model universe (here, primarily 
due to modified-gravity effects) has an effective equation-of-state parameter wx which is found to 
evolve toward the value —1 from above. 
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I. INTRODUCTION 



The fundamental theory of the strong interactions is nowadays taken to be quantum 
chromodynamics (QCD); see, e.g., Refs. |2| and other references therein. In the framework 
of this theory, there is evidence for the existence of a gluon condensate [s-G]. The question, 
then, is how the gluon condensate gravitates and evolves as the Universe expands. Here, a 
tentative answer is obtained by use of the so-called g-theory approach for the gravitational 
effects of vacuum energy density [z iO]. 



The outline of this article is as follows. In Sec. [TTl an example of a gluon-condensate- 
induced modification of gravity is presented and the corresponding field equations are de- 
rived, which are then reduced for the case of a spatially fiat Friedmann-Robertson- Walker 
universe. In Sec. IIIH the resulting evolution of a simple three-component model universe 
is studied both analytically and numerically, in order to establish whether or not a model 
universe can be obtained which resembles the observed "accelerating Universe" [l2|. In 
Sec. HVt concluding remarks are presented. 



II. QCD SCALE MODIFIED GRAVITY AND COSMOLOGY 



A. Theory: Action and field equations 



It has been argued [10| that, in a de-Sitter universe with Hubble constant H, a QCD- 
scale vacuum energy density pv ~ |-fr|AQ(jj3 could arise from infrared effects of the gluon 
propagator. Since the de-Sitter universe has Ricci curvature scalar \R\ ~ and the 
particular gluon condensate q has energy scale q ~ Aq^j-,, one is led to consider the following 
modified-gravity action {h = c = 1): 



S. 



eflf 



d^xy/^ Kf{R,q) + e{q) + CMW 



f = R+h = R + r]K-^ 



1 3/4 



(2.1a) 
(2.1b) 



with gravitational coupling constant K = [IGnG)^^ > 0, dimensionless coupling constant 
?7 > [standard general relativity has rj = 0], energy density e(g) of the gluon condensate 
q{x), and matter field ipi^) [later on, this single matter component will be generalized to 
matter components]. The precise definition of the gluon-condensate variable q{x) in the 
context of QCD has been given in Ref. 10|, to which the reader is referred for details. In 
the following, q is simply assumed to be nonzero and is, in fact, taken to be positive. The 
relation between the gravitational constant G and Newton's constant Gn ISi uM ^^^^ t)e 
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discussed in Sec. IIIIBI Throughout, the conventions of Ref. [l5| are used, in particular, 
those for the Riemann tensor and the metric signature ( — h ++)• 

The field equations from (12. ip are fourth order and it is worthwhile to switch to the scalar- 
tensor formulation which has field equations of second order. The equivalent Jordan-frame 
Brans-Dicke theory 
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S. 



(BD) 



off 



d X 
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has action 

K{<f)R-U{<l),q)]+e{q)+CMW 



U 



-(1/4) {vyK')\qf/yil-<P) 



(2.2a) 
(2.2b) 



in terms of a dimensionless scalar field restricted to values less than 1 [0 would be greater 
than 1 for the rj < case not considered here] . The dependence of potential (I2.2b|) allows 
for the so-called chameleon effect p^], which will be briefly discussed at the end of this 
subsection.^ The proof of the classical equivalence of the actions (12. ip and (12. 2p . for 77 7^ 
and q ^ 0, is not affected by the presence of the g-fleld in the function / of (]2.1bp . See, e.g., 
Refs. 22N24I for details of the proof, which is straightforward and need not be repeated here. 



Anyway, the classical equivalence of (12. ip and (12. 2 p can be verifled directly by eliminating 
(j) from (I2.2ap . using its fleld equation R = dU/dcj) with U{(f)) given by fl2.2bp . 

At this moment, two remarks may be helpful to place the theory considered in con- 
text. First, the rigorous microscopic derivation of the effective action (12. ip remains a major 
outstanding problem, because only a rough argument has been given in the appendix of 
Ref. p^], where rj was called / (see also Ref. [25| for a general discussion of the physics 
involved and 26| for a heuristic argument). Awaiting this derivation, the main motivation 
of (12. ip is that it naturally gives the correct order of magnitude for the present vacuum 
energy density (see Ref. 10| and Sec. II Vp . Just to be crystal clear: the term h in fl2.1bp 
is, at present, purely hypothetical and the aim of this article is to explore its cosmological 
consequences, leaving aside its theoretical derivation. 

Second, the effective action (12. ip is only considered to be valid on cosmological length 
scales and additional nonstandard terms in f{R, q) can be expected to be operative at 



smaller length scales, relevant to solar-system tests and laboratory experiments [22|, |23|. 
Purely phenomenologically, the h term in fl2.1bp could, for example, be replaced by an 
extended term 



r^K~^\qf^\R\"^/{\qf^ + CK^\R\) 



(2.3) 



^ See also Ref. [20| for chameleon-type effects in a different context and Ref. 2]J for recent analytic and 
numerical work on the scalar profiles from compact objects, extending the original analysis of Ref. 18 1. 
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with constants < 77 ^ |C| < 1. This term h^^t vanishes as at large enough 

curvatures and, for 77 ~ 10~^ and \(\ ~ 1, is consistent with the relevant bound in Ref. 23 1 
based on the Eot-Wash laboratory experiment {2^. 

Returning to the action fl2.2p . the field equations are obtained from the variational prin- 
ciple for variations 6g^i, of the metric g^^, variations 6(f) of the Brans-Dicke field 0, and 
variations SA of the microscopic field A responsible for q condensate (see, in particular. 



Refs. 



loll). Specifically, the field equations are 

1 „ „., 1 ^ „„\ 1 ~ „., 1 



with the covariant derivative V^, the invariant d'Alembertian □ = VV^, the energy- 
momentum tensor TjJ*/ of the matter field ip, the integration constant /i, and the effective 
energy densities 

~ de 

dU 

U^U-q — . (2.5b) 

Two comments are in order. First, the reason of having the extra term —qde/dq in f l2.5ap 
and —qdU/dq in ( I2.5bp is the fact that the field q is not fundamental but contains, in 



addition to the microscopic field A mentioned above, the inverse metric g^" (see Sec. II 



of Ref. [lOj). Second, the constant /i on the right-hand side of fl2.4cp can be interpreted, 
for spacetime-independent q and dU/dq = 0, as the chemical potential corresponding to the 
conserved charge q (see, in particular, the detailed discussion in Sees. II A and B of Ref. 

For completeness, also the generalized Klein-Gordon equation is given, which is obtained 
by taking the trace of fl2.4al) and using (]2.4bp : 



□ ^._L(r„-4T)+l&-i^|, (2.6) 

with the matter energy-momentum trace Tm = g^^. 
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Eliminating qdU/dq from f l2.4ap and fl2.4cp . the final field equations are 

R-''-^Ra-'' = -^{TM-Pv3-'')-j^Ucr''-^[\/'-\/''-gi"'ny, (2.7a) 

in terms of the gravitating vacuum energy density 

pv(g)=e(g)-/ig, (2.8) 

with the integration constant /i. Equally, the generalized Klein-Gordon equation (12 .Gp 
becomes 

□ ^=5i.(r„-4p.)4c/-i^^. (2.9) 

where the very last term on the right-hand side, in particular, is relevant to the previously 
mentioned chameleon effect. With ( I2.7bp . this last term of (12. 9p becomes {—R/3)(f) and 
corresponds to an effective mass square term for the scalar field, with a mass square of the 
order of pu/K for the case of a pressureless perfect fluid. This is indeed one aspect of the 
chameleon effect, namely, an effective mass value dependent on the environment [l^. 



B. Differential equations for a flat FRW universe 



For a spatially fiat {k = 0) Friedmann-Robertson-Walker (FRW) universe [15| with 
scale factor a(r) and matter described by a perfect fluid, the 00 and 11 components of the 
generalized Einstein field equation (I2.7ap can be combined to give a generalized Friedmann 
equation. Together with equations obtained directly from fl2.7bp and (12. 9p . the relevant 
equations are then 

H'cP = -^p,ot-]:U-H>P, (2.10a) 
b K D 

H = -2if^-^^, (2.10b) 

0(p 

= _3i70 + _L(^p^^^_3P^^^)_^f/ + i0^, (2.10c) 

with the overdot standing for the derivative with respect to r (the somewhat unusual no- 
tation r is used for the dimensionful cosmic time, in order to reserve the letter t for the 
dimensionless time later on). The total energy density and pressure are given by 

Ptot = Pv + PM, Ptot = -Py + -Pm , (2.11a) 
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for the gravitating vacuum energy density 



Pv(g) = -Py(g) = e(g)-/ig, (2.11b) 

as discussed in the previous subsection. Observe that fl2.10ap reproduces the standard 
Friedmann equation for f/ = 0, = 1, and K = (IGnG)^^ = (IGttGtv)^^ = Kj^. 

The last two equations in f l2.10p are, respectively, first- and second-order ordinary differ- 
ential equations (ODEs) for H and 0. Two further ODEs can be obtained as follows. First, 
multiplying f l2.7cp by q gives an equation for the time dependence of the vacuum energy 
density, 

Pv = K (^U-<P^^ , (2.12a) 

which describes the energy exchange between the vacuum and the nonstandard gravitational 
field {U 7^ 0). Second, the standard energy conservation of matter gives 



Pm 



-3H [pM + Pjv/) = -^H {l + wm) Pm , (2.12b) 



where the matter equation-of-state (EOS) parameter wm = Pm/pm has been introduced 
(henceforth, wm will be assumed to be time independent). Equation fl2.12bp implies that, 
for the theory considered, there is no energy exchange between vacuum and matter (such 
an energy exchange for a different version of g-theory has been studied in Ref. {28|). 



C. Dimensionless variables and ODEs 



Now rewrite the cosm olog ical equations in appropriate microscopic units. The gluon con- 
densate q from Refs. js], llo| has the dimension of energy density, [q] = [e], which implies 
that the corresponding integration constant p is dimensionless, [p] = [1]. The equilibrium 
value go of the gluon-condensate variable q is taken to be determined by a laboratory experi- 
ment in an environment with negligible spacetime curvature and has the order of magnitude 
go = Eqcd — 0(10^ eV^); see Sec. IIII CI for further remarks. From this moment on, consider 

matter components, labeled by an index n = 1, . . . , N. 

Specifically, the following dimensionless variables t, h, f, r, u, and s can be introduced: 
r ^ t K/ql" , H{r) ^ h{t) q^ / K , (2.13a) 

g(r) ^ f{t) go , p(r) = r{t) ql'^ / K , (2.13b) 

f/(r) = u{t) ql'^/K^ , 0(r) = s{t) . (2.13c) 
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Observe that all dimensionless quantities are denoted by lower-case Latin letters. A further 
rescaling t = t'/ri and h = h' rj will not be used in the present article, as the effects from the 
unknown coupling constant rj are preferred to be kept as explicit as possible. 

It is, then, straightforward to obtain the dimensionless versions of the algebraic equation 
(I2.7cp . the last two ODEs in fl2.10p . and the matter conservation equation ( I2.12b|) gener- 
alized to matter components. This gives a closed system of 4 + equations for the 
4 + A^ dimensionless variables f(t), h(t), s(t), v{t), and rM,nit). Specifically, this system of 
equations consists of a single algebraic equation, 

drvif) _ dujsj) 

df - df ^ ^'■'^> 

and 3 + A^ ODEs, 

h = -2h^ — , 2.15a 
6 OS 

s = V, (2.15b) 

1 / . , 2 1 du , , 

V = - [rtot-^Ptot) -^hv - -u + -s — , (2.15c) 
6 ^ ^ 3 3 as 

rM,n = -3h [l + WM^n) rM,n, (2.15d) 

where, now, the overdot stands for differentiation with respect to the dimensionless cosmic 
time t and the dimensionless total energy density and pressure are given by 

N 

not = +'rv + ^ rM,n, (2.16a) 

n=l 
N 

Ptot = -rv + ^ u;M,nrM,n, (2.16b) 

n=l 

with matter EOS parameters WM,n still to be specified. The dimensionless vacuum energy 
density ry appearing in the above equations will be discussed in Sec. IIIDI The dimensionless 
potential u has already been defined by f l2.2b|) and fl2.13cl) . but will be given again in Sec. Ill Dl 
With the solution of Eqs. f l2.14p -( l2.15p for appropriate boundary conditions, it is possible 
to verify a posteriori the Friedmann-type equation fl2.10ap in dimensionless form: 

h^s + hv = {rtot-u)/6, (2.17) 

which, in general, is guaranteed to hold by the contracted Bianchi identities and energy 



conservation (cf. Refs. |l5|, l28|). Specifically, if the solution of Eqs. fl2.14p - fl2.15p satisfies 



f l2.17p at one particular time, then f l2.17p is satisfied at all the times considered. The 
additional constraint 02.1 7p will provide a valuable check on the numerical solution of the 
equations. 
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D. Ansatz for rv{f) and solution for f{s) 



The only further input needed for the cosmological Eqs. f l2.14p -( l2.15p is an Ansatz for 
the gravitating vacuum energy density pv(q) from (12. 8p or the corresponding dimensionless 



quantity ry from (12.13bp . In Refs. j? 



10| , it was argued that the vacuum variable q of the late 
Universe is close to its flat-spacetime equilibrium value go and the quadratic approximation 
can be used 

ry =7(1-/)', (2.18) 

with positive constant 7. 

From the ry definition in (12.13bp . the constant 7 in fl2.18p can be expected to be of order 



Z ^, with definition 

Z = q'J' ~ IGtT (^QCD/^Planck)' ~ 10"=^^ , (2.19) 

for the quantum-chromodynamics energy scale -Eqcd ~ 0.2 GeV and the standard grav- 
itational energy scale -Epianck = ^/hd' /Gn ~ 1-22 x lO"*^^ GeV (having set G ~ G^] see 



Sec. IIIIBp . According to the discussion in Refs. [l lO], / can also be expected to be suf- 
ficiently close to 1, in order to reproduce an ry value of order unity or less for the present 
Universe. For technical reasons, the value Z = 10~^ is taken in a first numerical study 
(Sec. nil Cp . Later, the proper boundary conditions and scaling behavior are considered 
rSec. HTTDD . 

The dimensionless scalar potential u{s, /) from fl2.2bp and fl2.13cp can be written as 
u{t) ^ UK'q,"' = -(r^V4) (2.20) 



where a relatively small value for rj appears to be indicated [10| by the measured value of 
the vacuum energy density; see Sees. IIIIBI and IIIIDI for further discussion on the numerical 
value of 1]. 

With the specific functions fl2.18p and fl2.20p . Eq. fl2.14p is a quadratic in a// and the 
positive root gives 



f^{s) = i^,/l + D{sy±D{s)J , (2.21a) 
D{s) = - s| > 0, (2.21b) 

K = (3/32) 7^77 > 0, (2.21c) 

where the minus sign inside the outer parentheses on the right-hand side of fl2.21ap holds for 
s < 1 [the plus sign appears for the s > 1 case not considered here]. Expression (I2.21ap can 
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then be used to eliminate all occurrences of / in the 3 + ODEs fl2.15p for the remaining 
3 + variables h(t), s(t), v(t), and rM^nif)- Referring to the ODEs (12.151) in the following, 
it will be understood that / has been replaced by f _{s) from f l2.2ip . 

III. THREE-COMPONENT MODEL UNIVERSE 
A. Preliminaries 

The modified-gravity theory considered in this article has been presented in Sec. Ill Al and 
the corresponding dynamical equations for a spatially fiat FRW universe in Sees. Ill BHll Dl 
The specific model studied in this section is a simplified version with only three components 
labeled n = 0, 1, 2: 

0. A gluon condensate [described by the dimensionless variable /] with dimensionless 
energy density ry{f) from f l2.18p and constant equation-of-state parameter wy = — 1, 
which is taken to give rise to a nonanalytic term in the modified-gravity action (12. ip . 

1. A perfect fluid of ultrarelativistic matter [e.g., photons] with energy density r^/,! and 
constant EOS parameter wm,i = 1/3. 

2. A perfect fluid of nonrelativistic matter [e.g., cold dark matter (CDM) and baryons 
(B)] with energy density rM,2 and constant EOS parameter wm,2 = 0. 

From the scalar-tensor formalism of the gluon-condensate-induced modification of gravity, 
there is also the auxiliary Brans-Dicke scalar s{t) to consider, with the dimensionless po- 
tential u{s, f) from ^TMf . 

The relevant ODEs follow from (I2.15P by letting the matter label run over n = 1,2. 
The ideal starting point of the calculations would be some time after the QCD crossover 
at T ~ Aqcd with rM,i ^ '"m,2- The physical idea is that the expansion of the Universe 
was standard up till that time and that, then, a type of phase transition occurred with the 
creation of the gluon condensate. Clearly, the gluon condensate can be expected to start 
out in a nonequihbrium state, / 7^ 1 and s 7^ 1. These issues will be discussed further in 
Sec. lITTDl 

At this moment, it is useful to recall the basic equations of a standard fiat FRW uni- 



verse 
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29| with gravitational coupling constant G = Gn or K = Kj^. For two components. 



a pressureless material fluid labeled M and an unknown fluid labeled X, these equations are 
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Qh"^ = Q[(x/af = tm + rx , (3.1a) 

-12 a/a = tm + rx + ^Pm + ^Px = tm + rx (l + 3wx) , (3.1b) 

where pm in fl3.1bp has been set to zero and the EOS parameter wx = Px/^x has been 
introduced. The standard energy- density parameters are defined as follows: 

= rM/(6 h^) , nx = rx/iQ h'') = 1 - • (3.2a) 

In addition, the following combination of observables can be introduced to determine the 
unknown EOS parameter: 

_2/aal\l /or^i\ 



where the last equality holds, again, for pm = 0. See, e.g., Refs. [30|, l3l| for details on how 
to reconstruct the dark-energy equation of state from observations. 
In order to be specific, take the following fiducial values: 



{l^M, nx. W^}""^'" = {0.25, 0.75, -1} , (3.3) 

which agree more or less with the recent data compiled in Refs. 32h37||. The standard fiat 
FRW universe with parameters (13.31) corresponds, in fact, to the basic ACDM model 2^ 
with CDM energy density tm o: 1/a? (with constant EOS parameter wm = 0) and time- 
independent vacuum energy density Z = rx oc a° (with constant EOS parameter wx = — 1 
and I the dimensionless version of the cosmological constant A). 

Returning to the modified-gravity theory fl2.ip - fl2.2l) . the same observables fl and wx can 
be identified. Specifically, the generalized Friedmann equation f l2.17p gives 

nx + ^M = 1, (3.4a) 

= l^grav + ^^y, (3.4b) 

Ogi-av = 1 — s — s/h — u/{6h'^) , (3.4c) 

Qv = rv/iQh^) , (3.4d) 

Qm = rM/{6h^) , (3.4e) 

where figrav is the new ingredient, as it vanishes for the standard theory with u = and 

s = 1. Similarly, the effective EOS parameter of the unknown component X can be extracted 

from mm and (1^X7]) for pm = 0: 

__2/aa 1\ 1 ry — u — Ahs — 2 s 

^ "3 V(^^2j Y^^i^~ rv-u-&hs + rM{l-s)) ' ^^'^^ 
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The right-hand side of (13. 5p shows that wx of the modified-gravity model (12 ■2p approaches 
the value —1 in the limit of vanishing matter content and constant Brans-Dicke scalar s 
as t — > oo. A priori, there is no reason why this approach cannot be from below, so that 
1 + Wx would be negative for a while (cf. Ref. 38|). 

The main goal of this section is to get a quasirealistic model for the "present universe," 
which is taken to be defined by a value of approximately 0.25 for the matter energy-density 
parameter Qm- This can only be done with a numerical solution of the ODEs, but, first, 
analytic results relevant to the asymptotic behavior at early and late times are discussed. 



B. Analytic results 

It is not difficult to get two types of analytic solutions of the combined ODEs f l2.15p and 
( I2.17P for the specific functions fl2.18p and f l2.20p . having used solution fl2.2ip to eliminate 
/ in favor of s. The first corresponds to a Friedmann universe with relativistic matter and 
without vacuum energy. The second corresponds to a de-Sitter-type universe without matter 
and with an effective form of vacuum energy. 

For rj = 0, the first analytic solution of f l2.15p -f l2.2ip has only relativistic matter {wm,i = 
1/3) contributing to the expansion. Specifically, this Friedmann solution (labeled "F") is 
given by 

/i(^) = (1/2) r\ s(^) = /(F) = l, (3.6a) 

rZ\ = (3/2) t~' , rJ), = . (3.6b) 

Remark that standard general relativity [which has, from the start, the action equal to (12. ip 
for ?7 = and G = Gx] allows for arbitrary values rM,i(l) and rM,2(l) at reference time 
t = 1. 

For rj > 0, the second set of analytic solutions of fl2.15p -f l2.2ip has only vacuum energy 
contributing to the expansion, together with the effects of the gluon-condensate-induced 
modification of gravity {wx = — 1)- This type of solution has constant (time-independent) 
variables h > and s G (0, 1), with / given by f ]2.21ap . From f l2.15ap and f]2.15cp . using 
f l2.20p . a cubic in s is obtained, which needs to be discussed first. 

Specifically, the cubic in x = 1 — s reads 

9x^-6x2+ (1 + 9^2) x-6k2 = 0, (3.7) 

with parameter k defined by fl2.21cp . This cubic has three distinct real solutions for < k,^ < 
(5 v^-ll)/18 ^ (0.100094)^ Two of these solutions (with 2/3 < s < 1) give stationary de- 
Sitter-type solutions of the ODEs fl2.15p - fl2.2ip . These two roots can be written in manifestly 
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real form by use of the Chebyshev cube root 



\t\<2 



2 COS [(l/3)arccos(t/2)] 



Ci/3(0) = v^. 



(3.8a) 
(3.8b) 



(3.9a) 
(3.9b) 



Defining the auxiliary parameters 

p = (1/3) {1/27 + 

q = (2/9)(l/82-2«:2),, 
the relevant roots of (13.71) are 

Shigh = 7/9 + VP Ci/3(- qp-'/^) , (3.10a) 

Smid = 7/9 + ^ [Ci/3{qp-'^') - Ci/3(- qp-'^')] , (3.10b) 
where the third solution siow = 7/3 — Shigh — s^i^ can be omitted, as it lies below 2/3 for k, 



in the domain considered [the stationary limit of, e.g., Eq. fl2.15cp requires s > 2/3 because 
ry from f l2.18p is non-negative by definition]. 

The first de-Sitter-type solution (labeled "deS,0" because / ~ for ^ 1) is then given 

by 

e(deS,0) 



/ 



(deS,0) 



Shigh = 1 - 6 - 162 + O(k^) , 
7_(shigh) =9(/s:2 + 36k^) +0(ft:'') 



(3.11a) 
(3.11b) 



/i(des,0) ^ n/{AV3) \f 



(dcS,0)|3/4 



,{deS,0)| 



^7/6 



[l-(81/2) «:^ + 0(/€6)] 



JdeS,0) 



0, 



(3.11c) 

in terms of the function /_(s) defined by (12. 2 lap and with an integer n = 1,2 to label 
the different matter components. Note that the expression in the middle of fl3.11cp simply 
follows from fl2.15ap for h = and u from (12.200 . 

The second solution (labeled "deS,l" because / ~ 1 for ^ 1) is given by 



,(deS,l) 



Smid = 2/3 + k + 3k^ + (27/2) + 81 + O (k^) 



j(des,i) _ f_ (3^^.^) = 1 _ 6 K - 27 - 162 /t^ + O , 

/,(deS,l) ^ ^/(4v^) |/(deS,l)|3/4 | ^ _ ^(dcS.l) | -1 ^ ^/2^/l024 



1024 - 1536 K + 1152 /t^ + 1728 + 17496 + 0{k^) 



„(deS,l) 



(3.12a) 
(3.12b) 

(3.12c) 
(3.12d) 
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where k is non- negative according to the original definition fl2.21cp . Note that the last 
expressions of both fl3.11cp and fl3.12cp are proportional to with all further dependence 
on 7 entering through the parameter k oc ^7^/7, as can be expected on general grounds from 
the ODEs flZT^ without matter. 

It is not quite trivial that there indeed exist de-Sitter-type solutions in the modified- 
gravity theory (12. ip . The first solution (13. lip is far from the equilibrium state /equii = 1 and 
the second solution (I3.12p is close to it, at least for \n\ <^ 1. The scaling behavior of both 
solutions under the limit 7 —> 00 for constant rj is also different, with h diverging for the 
first solution and staying constant for the second. For fixed parameters 7 and 1], numerical 
results suggest that the first solution (I3.1ip is unstable and the second solution (I3.12p stable 
[and possibly an attractor]. In the following, the focus is on the second solution close to the 
equilibrium value feqmi = 1 (corresponding to g = go)- 

In fact, two remarks on the de-Sitter-type solution (I3.12p are in order. First, observe that 
local experiments in this model universe with (p^^'^^'^^ ~ 2/3 < 1 would have an increased 
effective gravitational coupling 

(dcS,l) 



^ local cxps 



(l/(f)id-s,i)^G^{3/2)G, (3.13) 



where the term G / (p^'^^^'^^ in the middle comes directly from the combination K (f) = 
(p/llGnG) present in the action (12. 2p . Here, "local experiments" denote experiments on 
length scales very much less than the typical length scale of de-Sitter-type spacetime, the 
horizon distance Lhor = cH^^'^^'^\ whose numerical value will be discussed shortly. It would 
then appear that the quantity (I3.13P must be identified with Newton's gravitational con- 
stant Gtv as measured by Cavendish [3| and modern-day experimentalists [l^; see 39] for 
additional comments. 

Second, the de-Sitter-type solution (I3.12p of model (12. 2p or equivalently model (12. ip has 
the inverse Hubble constant 

(;,(des,i))-i = 4/^3 ^-1 ^ 2.3 X 10^ (^^^ , (3.14) 

as follows from (I3.12cp by neglecting terms suppressed by powers of k = 0(1/7) = O(10~^^) 
and anticipating a particular order of magnitude for the model parameter 77. With the 
conversion factor from (I2.13ap . the dimensionless quantity (I3.14p corresponds to 

~ 4/73 (3/2) A-„ ~ 8 X 10" s (15^) [^^J . (3.15) 

where, according to (I3.13p . an approximate factor 3/2 appears in going from K to the 
Newtonian value = {IGttGn)^^ . The time scale found in (I3.15p is of the same order as 
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the inverse Hubble constant (Hri) ^ ~ 4.5 x 10^^ s {0.70/ ho) for the measured value Hq ^ 0.70 



as reported in Refs. [32|, |36|, |37 |. 



By equating the theoretical quantity l/H^'^^^'^^ from (13.151) multiplied by an ad hoc factor 
g = with the measured value 1/Hq, a first estimate of the model parameter rj in the original 
action (12. ip is obtained, 



V 



VS Kj^ q-^^^ Ho 10-\ (3.16) 



for the go and Hq values mentioned in the previous paragraph. Admittedly, the choice of 
one-half for the factor g is somewhat arbitrary, but consistent with the physical picture of 
our present Universe entering a de-Sitter phase. A more reliable estimate of t] will come 
from the numerical study of a model universe with both vacuum and matter energies. The 
numerical solution found will be seen to interpolate between the analytic solutions (13. 6 p and 

(EIS]). 



C. Exploratory numerical results 

Equation (I2.15ap for the potential u{s, /) from (12.201) makes clear that a model uni- 
verse with an asymptotically nonvanishing Hubble constant, h{t) — )■ const 7^ 0, requires 
a nonvanishing modified-gravity parameter, 77 7^ 0. The analytic de-Sitter solution with 
h = s = f = has already been given in Sec. IIII B[ 

The numerical solution of ODEs (12.151) for rj ~ 10^^ is presented in Fig. [T] and several 
observations can be made: 

(i) The boundary conditions on the functions will be discussed in Sec. IIII D[ 

(ii) There is a transition from deceleration in the early universe to acceleration in the late 
universe. 

(iii) The values for s, 1 — /, and h at the largest time shown in Fig. [T] agree already at the 
10% level with those of the analytic de-Sitter- type solution (13.120 . 

(iv) The ratio rM,tot/ (6 h'^) is equal to 0.25 at the dimensionless cosmic time t !^ 1.4 x 10'^. 

Points (ii)-(iv) suggest that, for the model parameter values chosen, the model universe at 
tp = 1.432 X 10^ resembles our own present Universe, characterized by the values (13. 3p . 

More quantitatively, the following three estimates can be obtained. First, the product 
of the dimensionless age tp of the present universe with its dimensionless expansion rate 
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FIG. 1: Numerical solution of ODEs (j2.15p . with vacuum energy density (j2.18p . Brans-Dicke 
scalar potential (|2.20p . and both relativistic matter (energy density rjv/,i) and nonrelativistic 
matter (energy density rM,2)- The figure panels are organized as follows: the panels of the 
first column from the left concern the expansion factor a(t), those of the second column the 
modified-gravity scalar s(t), those of the third column the gluon-condensate vacuum variable /(t), 
and those of the fourth column the matter energy densities rM,n- The model parameters are 
(7, "rf, wm,1: wm,2) = (10^) 9 X 10~^, 1/3, O), with the resulting parameter k = (3/32) 77^/7 = 
8.4375 X 10^-*^^. The boundary conditions at tstart = 0.1 are (a, /i, s, f , 1 — /, rM,ii rM,2) = 
(1, 4.082483, 0.8, 0.8164966, 8.437500 x 10"^ 75.97469, 24.02531); see Sec. [mD] for details. The 
several energy-density parameters 17 and the effective "dark-energy" equation-of-state parameter 
wx are defined in (13. 4p and (13. 5j) . respectively. With 7//?^ ^ 1, the values of are negligible 
compared to those of ^grsv for the time interval shown. 

h{tp) ^ 0.6351 X 10-3 gives 

/i(tp) ^ 0.91 , (3.17a) 

which also holds for the product of the dimensionful quantities, TpH{Tp) ^ 0.91. 

Second, evaluating the particular combination (13. 5p of first and second derivatives of a(t) 
and the matter energy density pm, the present effective EOS parameter of the unknown 
component is found to be 



_ 2 / da 1\ 1 

Wx{tp) = -- + 



-0.66. (3.17b) 



3 V(a)2 2 J l-^M 
For larger times t^tp, this parameter wx{t) drops to the value —1, as can be expected from 
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the right-hand side of (13. 5p . Additional numerical values are wx = —0.75082, —0.98921, 
-0.99780, and -0.99989 for t = 2000, 4000, 8000, and 16000, respectively. Observe that 
the particular combination of observables (13. 5p is designed to be interpreted as the effective 
EOS parameter of the unknown component X only if matter-pressure effects are negligible 
(t > 500 in Fig. [ID . 

Third, consider the transition of deceleration to acceleration mentioned in point (ii) above. 
In mathematical terms, this time corresponds to the nonstationary inflection point of the 
function a{t), that is, the value tmficct at which the second derivative of a{t) vanishes but 
not the first derivative. Referring to the model universe at tp = 1.432 x 10'^, the inflection 
point tinfiect ~ 0.863 x 10^ corresponds to a redshift 

^inflect = a(ip)/a(^inflect) - 1 ~ 0.5 , (3.17c) 

which implies that the acceleration is a relatively recent phenomenon in this model universe. 
Inspection of the lower panels of Fig. [1] shows that the acceleration sets in when the ratio 
of Qx = ^grav + and ^jv/.tot is approximately unity, whereas the standard ACDM model 
would have fix/^A/,tot ~ 1/2 according to (13.1bp . 

Returning to the first estimate (I3.17ap . note that this quantity can be interpreted as the 
age of the present universe in time units obtained from the present expansion rate. But 
it is also possible to obtain the absolute age of the model universe, using the time scale 
contained in (I2.13ap . which requires as input the experimental value of the QCD gluon 
condensate go and the one of Newton's constant Gn, taken to be equal to the effective 
gravitational coupling Gn from (I3.13p . With the conversion factors from (I2.13ap and the 
relation G ~ s(tp)GN for K = l/(167rG), the numerical results tp ~ 1432, h{tp) ^ 1/1575, 
and s{tp) ~ 0.7267 give the following two dimensionful quantities of the present universe: 

Tp = tpKq~^'^ ~ 13.1 Gyr, (3.18a) 

Hp = h{tp) K-^ ql'^ ~ 68 km s"^ Mpc"^ , (3.18b) 

where the numerical values have been calculated with go = (210 MeV)^. Remark that, if 
the relation G ^ Gn holds for Cavendish-type experiments as mentioned in 39|, the same 



numerical values are obtained in (I3.18P by taking go ~ (190 MeV)'* and, if G ~ Gn /'^ holds, 
by taking go ~ (230 MeV)"^. All of these three go values lie below the value go ~ (330 MeV}^ 
indicated by particle physics j^], but the uncertainty in the latter value appears to be large {4- 



. In addition, it may be that certain particle-physics experiments are more appropriate 
than others to determine the truly homogeneous condensate go relevant to cosmology. 
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Next to the observations 



0.Q 



32 



37] . the values obtained in fl3.17p and fl3.18p have 



the correct order of magnitude, which is all that can be hoped for at the present stage. Still, 
it is remarkable that more or less reasonable values appear at all 4ol |. 

For comparison, the standard fiat-ACDM model fl3.ip - fl3.3p with boundary condition 
rM(tp)/rv = 1/3 gives the product TpH{Tp) ^ 1.01, the effective EOS parameter wx = — 1, 
and the inflection-point redshift 2;inflect = (6)^/'^ — 1 ^ 0.82. These three numbers fit the 
observational data perfectly well, but the ACDM model is purely phenomenological and 
cannot explain, without further input, ^ the absolute age of the Universe as in (I3.18ap or the 
absolute vacuum energy density as will be discussed in Sec. HVl 



D. Elementary scaling analysis 

In the previous subsection, the ODEs (12.150 have been solved numerically for certain 
parameter values and boundary conditions at t = tstart? which need to be discussed further. 
As explained in Sec. IIII At tstart is considered to correspond to a time just after the QCD 
crossover has happened. This implies, in particular, that the starting value /i(tstart) for the 
expansion rate is approximately given by the value [(ry + rA/.tot)/6]~^''^ of the corresponding 
standard FRW universe (13. lap . The / value at tstart follows from (12.2 ip for the chosen s 
value (see below) and the starting value for v is obtained by solving (12.170 . considered as a 
linear equation in v with all other quantities given. 

Next, the value of tstart itself and the corresponding values for r^,! and rM,2 need to be 
specified. These values depend on the physical ratio Z defined by (12.190 . Following the 
results for the standard FRW universe, take 

7 = 7 (3.19a) 

tstart =t^/Z, (3.19b) 

rAf,i(Wt) = rZ-'/{l + , (3.19c) 

rM,2(Wt) = rZ-'/y{l + Z'/') , (3.19d) 

where the constants 7, t, and r are numbers of order unity [in the present elementary analysis, 
they are just set equal to 1]. With t = 1 and the particular Ansdtze (I3.19cp - fl3.19dp . there 



^ Taking as additional input the measured value 32| ho « 0.70 of the Hubble constant Hq 



ho 100 km s~i Mpc ^ = ho (9.778 x 10^yr)"\ the ACDM-model result tq i?o « 1-01 gives the dynamic 
age To « 14.2 Gyr. 
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TABLE I: Numerical results for the "present epoch" [defined by Ojv/ (tp) = 0.25] in model universes 
with different numerical values for the parameters Z and r/, where the latter parameter controls the 
modified-gravity term in the action (|2.1|) and the former is defined by (|2.19p in terms of the physical 
energy scales. Other parameters and boundary conditions are given by (j3.19p . with constants 7, 
t, and f set equal to 1. A further boundary condition is s(tstart) = 0.8 ; see Sec. IIIIDI for details. 
The effective equation-of-state parameter wx and the inflection-point redshift ^inflect are defined 
in ()3.17bp and (j3.17cp . respectively. Figure [T] for Z = 10^^ illustrates the general behavior of h{t), 
wxit), and other physical quantities. 



z 


106 ^2 


10-3 tp 


10"^ h{tp) 


s{tp) 


tp h(tp) 


Wx{tp) 


•^inflect 


10-1 


0.8 


1.522 


5.980 


0.7272 


0.910 


-0.669 


0.541 


10-2 


0.9 


1.432 


6.351 


0.7267 


0.910 


-0.662 


0.538 


10-4 


0.7 


1.629 


5.584 


0.7259 


0.910 


-0.663 


0.515 


10-8 


0.8 


1.523 


5.967 


0.7255 


0.909 


-0.660 


0.505 


10-16 


0.9 


1.436 


6.330 


0.7256 


0.909 


-0.660 


0.506 



is equality of the relativistic (label n = 1) and nonrelativistic (label n = 2) energy densities 
around t ~ 1, which is not entirely unrealistic if the present universe has t ~ 10^. 

Finally, the boundary condition value s(tstart) is taken between and 1. The results are, 
however, rather insensitive to the precise value of s(tstart); see 43] for selected numerical 
results. The explanation is that, independent of the precise starting value, s{t) increases 
rapidly until, at t ~ 1, it bounces back from the s = 1 "wall" and, then, slowly descends 
towards the de-Sitter value, with some initial oscillations. 

Having specified the boundary conditions of the physical variables, the optimal model 
parameter 77 needs to be determined. The strategy is as follows: for a given Z value, assume 
an Tj value, determine tp from the condition fiM,tot(^p) = 0.25, evaluate the product tph{tp), 
and, if necessary, return to a new value of rj in order to get tp h{tp) closer to the asymptotic 
value of approximately 0.909. 

Numerical results are given in Table HI Three physical quantities, the relative age of the 
present universe tph{tp), the effective EOS parameter Wx, and the inflection-point redshift 
-2^inflect, appear to approach constant values as Z drops to zero. This nontrivial result suggests 
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that the behavior shown in Fig. [T] and the corresponding estimates (13.17p -( l3.18p also apply 
to the physical case with Z ~ 10^^^ as given by f l2.19p . 

IV. CONCLUSION 

The bottom-row panels of Fig. [H if at all relevant to our Universe, suggest that the 
present accelerated expansion may be due primarily to the nonanalytic modified-gravity 
term in the action fl2.ip rather than the direct vacuum energy density pvil), because q is 
already very close to its equilibrium value go, making pviq) ~ Pvi.%) = 0. Referring to the 
definitions in fl3.4p . the second panel of the bottom row shows the effective energy- density 
parameter f2grav due to the gluon-condensate-induced modification of gravity and the third 
panel the energy-density parameter VLy from the vacuum energy density proper [with EOS 
parameter wy = —1], their total giving VLx which equals I — VLm for a flat FRW universe. As 
discussed in Sees. IIII Al and IIII C[ the total unknown 'X' component has an effective EOS 
parameter wx which drops to the value —1 as the de-Sitter -typ e universe is approached. 



Remark that, in contrast to the results of, e.g., Refs. [22|, |23|, nontrivial dark-energy 
dynamics has been obtained, because the effective action f l2.ip is assumed to be valid only 
on cosmological length scales, not solar-system or laboratory length scales [see also the 
discussion in the paragraph of Sec. Ill Al containing Eq. f l2.3p ]. As it stands, the effective 
action (12. ip can be viewed as an efficient way to describe the main aspects of the late 
evolution of the Universe, with only two fundamental energy scales, -Eqcd ~ 10^ eV and 
-Epianck ~ 10^*^ eV, and a single dimensionless coupling constant, 7] ~ 10^'^. Moreover, this 
effective coupling constant rj can, in principle, be calculated from quantum chromo dynamics 
and general relativity, which may or may not confirm our numerical value of approximately 
10-3; cf. Refs. 
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25| and the third remark in the Note Added. 
Elaborating on the source of the present acceleration, consider the second term on the 
right-hand side of (I2.7ap . which can be rewritten as +{2(j)K)~'^ (pv,bd) Q^v for the Brans- 
Dicke vacuum energy density Pv,bd = —KU. The exact de-Sitter-type solution (I3.12p for 
K ^ 1, together with the conversion factor from (I2.13cp and Newton's constant from (I3.13p . 
then allows for the following estimate: 



PV,BD 



(deS.l) 



-u qf/K • ■ = 127r r,' q^ G ~ (vr/S) r,' i^^cD/^Pianck 

~ (2 X 10-3 eV)^ X i-^X { , (4.1) 
^ ' VIO-V l^(420MeV)V ^ ' 

where go has been expressed in terms of the QCD string tension -R'qcd fl, specifically, 



(dcS.l) 



.2 ^3/2 
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go = -E-QCD ~ The parametric dependence of the above expression, pv (x 

-^QCD /-^Planck) the Same as that of the previous estimate (6.7) in Ref. jlo|, but expression 
(14. ip now comes from the solution of field equations. Two other dimensionful quantities, 
the age and expansion rate of the Universe, have already been given in fl3.18p . 

Before the asymptotic de-Sitter-type universe with effective energy density (14. ip is 
reached, the Brans-Dicke scalar evolves and allows for an effective EOS parameter wx 
different from —1 [the scalar has no direct kinetic term in the action (I2.2ap . but the (pR 
term does give, by partial integration, an effective kinetic term for 0, which, in fact, leads to 
the generalized Klein-Gordon equation (12. 9p ]. For the present Universe, the general lesson 
may be that the deformation of the QCD gluon condensate q by the spacetime curvature 
of the expanding Universe can result in an effective EOS parameter wx which evolves with 
time and, for the present epoch, can still be somewhat above its asymptotic value of —1. In 
turn, a possible discovery of a wx time dependence may provide an additional incentive to 
theoretical investigations of the physics of the gravitating gluon condensate. 
Note Added. — After completion of the work reported here, we became aware of two earlier 
articles and a third article recently posted on the archive. The first article is a systematic 
study of the cosmology of f{R) modified-gravity models and identifies the modified-gravity 
term (I2.1bp . for constant q, as cosmologically viable [observe the different sign definition of 
R compared to ours]. The second article [45| investigates the growth of density perturba- 
tions in f{R) modified-gravity models and establishes, in Eq. (42), the effective gravitational 
coupling parameter for subhorizon CDM density perturbations, which turns out to be close 
to Gm for the model universe of Fig. [T]at times t < 500 (redshifts z > 1). The third arti- 
cle 46| presents a QCD calculation for the origin of the modified-gravity term (]2.1bp and 



may also explain the smallness of the coupling constant rj, even though many conceptual 
and technical issues remain to be resolved. 
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